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Abstract 

An ID tight-binding version of the Dirac equation is considered; 
after checking that it recovers the usual discrete Schrodinger equation 
in the nonrelativistic limit, it is found that for two- valued Bernoulli po- 
tentials the zero mass case presents absence of dynamical localization 
for specific values of the energy, albeit it has no continuous spectrum. 
For the other energy values (again excluding some very specific ones) 
the Bernoulli Dirac system is localized, independently of the mass. 



In one-dimensional quantum systems general random potentials induce 
localization and no conductance PjJ |2] , irrespective of the disorder intensity. 
Exceptions are restricted to random models with local correlations, as poly- 
mer models random palindrome models @J, both including the important 
precursory random dimers E] (see |3] for rigorous approaches). In this 
article it is presented a random discrete model with no local correlation for 
which delocalization occurs in some situations; to the best of authors' knowl- 
edge it is the first tight-binding model with such property (in jH] dynamical 
delocalization is shown for a tight-binding Schrodinger model, but the ran- 
dom potential is decaying). Since mathematical proofs will appear elsewhere, 
it is hoped that the present note will contact physicists with recent interesting 
mathematical results on delocalization in one- dimension. 

The model is in fact very simple. It is a relativistic version of the well- 
known tight-binding Schrodinger Hamiltonian (with h — 1) 

(H S l/)) n = -77— (A^/On + V n 1p n = — (-1p n+1 - 4>n-l + tyn) + Vn^n (1) 
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(recall that it is common to take +A instead of —A, and also exclude the 
constant factor "2", in the kinetic term). For very general random potentials 
V n the model is localized, including the Bernoulli potential for which 
the site energy V n is assigned one of two values ±v at random (— v with 
probability < p < 1 and +v with probability 1 — p, say); the spectrum of 
the corresponding operator has no continuous component [2]. 

By imposing the strong local correlation that the site energies V n are 
assigned for pairs of lattices, i.e, Vm = V^+i = ±f for all n, one gets the 
random dimer model exhibiting delocalized states [3 EE IS]- Despite of 
dynamical delocalization the dimer Schrodinger operator has no continuous 
component in its spectrum. The first example of system with dynamical 
delocalization and pure point spectrum was a peculiar almost periodic op- 
erator jHj. Such results clarified the difference between mathematical local- 
ization (i.e., pure point spectrum) and dynamical localization (i.e., bounded 
moments, see ahead). In the zero- mass case, the Dirac [TB] model discussed 
here has pure point spectrum and dynamical delocalization with no added 
correlation to the Bernoulli potential. A crucial first step for the arguments 
will be the appropriate way of writing the transfer matrices and their simi- 
larity with those of the Schrodinger dimer model. 

Previous works have considered the one- dimensional Dirac equation and 
relativistic effects on conduction in disordered systems ^U], localization ^T] 
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and comparative studies of relativistic and nonrelativistic Kronig-Penney 
models with ^-function potentials [TTJ [T2|. Although tight-binding equa- 
tions for the electronic amplitudes have naturally arisen in some of these 
studies, the phenomenon reported here have not. The interesting question of 
the comparison of the relativistic and nonrelativistic localization length was 
numerically investigated ^T] in some cases, and it was found that which one 
is larger depends on the on-site energy and also on the energy particle. 

Consider a particle of mass m > in the one- dimensional lattice TL under 
the site potential V = (V n ); the Dirac tight-binding version is proposed as 

H D (m,c) = H (m,c) + VI 2 =(^ d ^ +mc 2 a 3 + VI 2 , (2) 

with c > being the speed of light, 03 the usual Pauli matrix, J 2 the 2x2 
identity matrix and d a finite difference operator (a discrete analogue of the 
first derivative) 

(#)„ = Vn+l - 4>n- 

Since d is not Hermitian, its adjoint (d*ip) n = ipn-i — 4>n appears in the 
definition of Hd (the inclusion of the imaginary unit % in front of the difference 
operators d and d* is immaterial). In case V n takes a finite number of values, 
it is clear that H & is a bounded Hermitian operator acting on £ 2 (2; (D 2 ) and 
the resulting Dirac equation can be recast in the compact form 

f„ ( W ,\ fmc 2 + V n cd* \ , , 

with the spinor * = (* n ) and = 

One can easily verify that the nonrelativistic limit of (j2J) is the equation 
associated to the Schrodinger operator (JTJ; this is an important support for 
the Dirac model just introduced. Following the traditional prescription for 
the nonrelativistic limit of the Dirac equation jTSj, first one removes the rest 

energy by inserting \1> = e ~ imc2 *$ = e ~ imc2t (®_ \ into © so that 



1— = c , , - 2mc z , _ + V®. 



dt \d 

For large values of c, the equation in the second row above can be solved 
approximately as <ft~ = d<p + /2mc, and inserting this into the first equation 
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results in 

^ = _W + ^. (4) 
Similarly, by considering \l/ = e tmc *<& one finds 

■ <A) 1 dd*<j)~ + v<jr. (5) 



9t 2m 

Since d*d = dd* = —A, then (jJJ and (0) correspond to the one-dimensional 
tight-binding Schrodinger equation associated to ((TJ with positive and nega- 
tive free energies, respectively. 

Another point directly related to the continuous Dirac equation is the 
presence of the so-called zitterbewegung jTH] phenomenon for (jHJ); here it will 
be explicitly considered the particular case of free particle and small mass m. 
Following Section 69 of Dirac's book [14J, let h denote the position operator 
(n 1 i r ) n = nty n , so that its time evolution under the free operator H (m, c) is 
n(t) = e lH Q t tie~ %Uat ] the velocity operator is then 

dn ( t ) _ .Tit A MM _ JH o t nA( n.\ e -iH t 



dt 



i[H Q ,h(t)\ = e lHot cA(0)e- lHot = cA(t), 



/ o —d* — 1 \ 

with A = A(0) = i f ^ ^ q J . Notice that A is Hermitian, A 2 = / 2 , 

so that its spectrum is ±1 and then the spectrum of cA is ±c. Since e~ lH ° l 
is unitary, it follows that the spectrum of the above velocity operator is ±c 
for all t. Hence it indicates that the possible speed measurements would 
result only in ±c. Now, for small mass m the time derivative of the velocity 
operator is given by 

= i[H , cA{t)\ = 2iH F(t), (6) 

with 

F(t) - l —u-^ p ^t [dd* \ iHot 

The operator F = F(0) anticommutes with H ; thus dF(t)/dt = 2iH F(t) 
and it is found that F(t) = e 2tHot F, which is fast oscillating. Inserting this 
into © one finds d(cA(t))/dt = dF(t)/dt; after integrating from to t one 
gets dh{t)/dt = cA — F + e 2lHot F and the velocity quickly oscillates around 
an average value; this is a version of zitterbewegung. 
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Now the localization results will be discussed. Consider model (J2J) with 
V n taking the values ±v , v > randomly. Denote by 5^ the elements of the 

canonical position basis of £ 2 (7L; (D 2 ), for which all entries are except at 
the nth entry, which is given by ( ^ J and (-^) f° r the superscript indices + 
and — , respectively. If \l/ n = " is a solution of the eigenvalue equation 



then it is simple to check that 



m 2 c 4 -(E-V n ) 2 
c 2 

mc 2 -E+V n 




is the transfer matrix at the nth step. Recall that the Lyapunov expo- 
nent "f(E) represents the average rate of exponential growing of the norm of 
transfer matrices 

\\rpE rpErpEU ^ 7(B) O. 

\\ 1 v„ 1 v 2 1 Vi II ~ e ) 

l/ , ~f(E) is called the localization length. A vanishing Lyapunov exponent is 
an indication of derealization, so the next task is to find possible energies E 
with j(E) = 0. In order to get vanishing Lyapunov exponents and diffusion, 
the arguments will follow those in jHJ [7j; detailed mathematical proofs will 
appear elsewhere [To] . 

Given an initial spinor \1/ with only one nonzero component (i.e., well- 
localized in space), the dynamical derealization will be probed by the time 
average of the mean squared displacement (also called second dynamical 
moment) 

M™{t) = - [Vn 2 (| (S+, e -tHn(m,c)s^ .2 + . / j- £ -iH D (m,c)s^ |2\ ds] 
t Jo ~ V ' 

dynamical localization is characterized by a bounded M^(t) < cte, for all t; 
otherwise the system is said to present dynamical derealization. 

First the spectral questions will be faced. By adapting the multiscale 
analysis HE] to the Dirac operator (J2J it is possible to show that, due 
to the random character of the Bernoulli potential, for typical realizations 
the spectrum of ifo(m, c) has no continuous component for any m > 0, and 
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with exponentially localized eigenfunctions. In other words, mathematical 
localization holds for Hp. 

By using as the main tool Furstenberg theorem , for m > it is shown 
that j(E) = if, and only if, E — (for v = c\J1 + m 2 c 2 ) and the four 
possibilities E = ±c/\^2 ± ca/2 + m?c 2 (for v = c/a/2); so, for other values 
of energies dynamical localization can be shown. For such E values with 
'j(E) = 0, it was not possible to give an answer about dynamical localization 
yet. 

Nevertheless, restricted to the massless (m = 0) case, if < v < c the 
Lyapunov exponent vanishes for E = ±t> and v ^ c/a/2, and following [3| it 
is possible to show that 

Ml{t) > ctet 3/2 , 

i.e., there is no dynamical localization despite the absence of a continuous 
component in the spectrum of the random operator Hd(0,c). Due to its 
importance here, it is worth including the main argument for the vanishing 
of j(E = v) (the case E = —v being similar). In this case, the possible 
transfer matrices are 




Notice that T^ v and are commuting matrices and both have spectral 
radius equal to 1 (for such v). If n_ denotes the average number of times 
that the potential — v occurs in n trials, then n — n_ is the average number 
of times that the transfer matrix is the identity. Thus, if p is the probability 
for the potential value —v, 

j(v)= lim -ln||T# ■••r^ll = lim — lnlKT^r-H 1 /™- =plnl = 0. 

The heuristics for the dynamical derealization in this case can be found in 
the paper by Dunlap, Wu and Phillips in Ref. 0. The main concern for the 
proof is the uniform boundedness of the product of transfer matrices [3 CH] • 
It is very important to stress that here derealization is not synonymous of 
zero Lyapunov exponent, as some people have considered. 

For small but nonzero mass, it is expected that the dynamics follow closely 
the massless case, at least for a small period of time. The final result to be 
reported is an inequality confirming such expectative; by making using of 
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Duhamel formula, it can be shown that, given an initial there exists 
C > so that, for all t > 0, 

<CmcH\ (7) 

Therefore, if the time t is not too large and/or the mass m is sufficiently 
small, the mean squared displacement follows rather closely the delocalized 
massless case, so that inattentive numerical simulations could give a wrong 
insight. 

It is natural that this model would be applied to any case the one- 
dimensional tight-binding Schrodinger operator was used; it would be the 
first step for their relativistic versions. Since the Dirac operator in the mass- 
less case presents dynamical derealization, this becomes a potential source 
for explaining some observed effects (at least for small m) as details in the 
theory of mesoscopic systems [TH] . 

Summing up, a natural one-dimensional Dirac tight-binding model was 
proposed which was supported by its nonrelativistic limit (it recovers the 
discrete Schrodinger model) and the presence of zitterbewegung. Then re- 
sults about mathematical and dynamical localization were reported for such 
operator with random Bernoulli potentials: for all values of c > and mass 
m > 0, there is mathematical localization, but in the massless case and po- 
tential intensity v < c, particular values of the energy imply the absence of 
dynamical localization, although no potential correlation was imposed. It 
is possible that this model is the simplest one with such derealization. Fi- 
nally, relation (J7|) gives quantitatively an estimate of how, for small time t, 
the dynamics of the localized regime follows the delocalized one. 
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